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Abstract. Density (or state) dependent pair potentials arise naturally from 
coarse-graining procedures in many areas of condensed matter science. However, 
correctly using them to calculate physical properties of interest is subtle and 
cannot be uncoupled from the route by which they were derived. Furthermore, 
there is usually no unique way to coarse-grain to an effective pair potential. Even 
for simple systems like liquid Argon, the pair potential that correctly reproduces 
the pair structure will not generate the right virial pressure. Ignoring these 
issues in naive applications of density dependent pair potentials can lead to 
an apparent dependence of thermodynamic properties on the ensemble within 
which they are calculated, as well as other inconsistencies. These concepts are 
illustrated by several pedagogical examples, including: effective pair potentials 
for systems with many-body interactions, and the mapping of charged (Debye- 
Hiickcl) and uncharged (Asakura-Oosawa) two-component systems onto effective 
one-component ones. 



PACS numbers: 61.20.Gy 



1. Introduction 

No known materials exist in nature whose behaviour can be completely captured by 
pair potentials alone. Even a noble gas like Argon has a finite contribution from 
three-body Axilrod- Teller triple-dipole interactions [Q. Thus, the pair-potentials used 
to describe condensed matter systems always arise from coarse-graining procedures, 
where a subset of the degrees of freedom of the full (quantum) statistical mechanical 
system are integrated out. In the aforementioned example, integrating the three-body 
interactions over angular coordinates results in effective parameters for the Lennard 
Jones (LJ) potential, which will depend on state. In metallic systems, integrating out 
the free electrons leads to a configuration independent volume term and pair potentials 
that depend on the global density 0. Alternate coarse-graining procedures for metals 
such as the embedded atom method[[| f|, effective medium theory[[|, Finnis-Sinclair 
potentials^] or glue potentials^ result in a local density dependence. 

Coarse-graining methods are crucial to deriving tractable statistical mechanical 
treatments of soft-matters systems, where a large number of different length and 
time-scales may coexist. An increasingly popular coarse-graining technique consists of 
deriving effective potentials and exploiting their analogy to well studied simple atomic 
or molecular systems to calculate phase behaviour and correlations ||, |9[ [n|. Again, 
these effective interactions are often reduced to an approximate pairwise description 
with parameters that depend on state. 

Direct inversion from experimental structure factors are another way to derive 
the parameters for effective pairwise potentials (O, 12j. These almost always show 
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a dependence on state, especially for the case of soft-matter systems. This is not 
surprising, of course, since one can easily imagine that the interactions between 
two effective particles depends on the overall density. For example, changing the 
concentration of a micellar solution may affect the internal structure of the micelles, 
which in turn leads to a density dependence of the effective pair interaction between 
the particles. 

That an effective pair potential derived in one context doesn't always perform well 
in another is well known, and usually categorised as a problem of transferability. For 
example, if the parameters of an effective pair potential depend on density (but the 
explicit density dependence is not known), then a parameterisation of the potential 
at pi is not the same as the one needed at a different density P2 - the potential at pi 
is not transferable to the state point at p%. 

What 1 will endeavour to show in the present paper is that there are deeper 
problems associated with the use of effective pair potentials, even when the problem of 
transferability appears to be solved. The particular example studied is pair potentials 
that depend on the global density p as v(r; The focus is partially pedagogical. For 
that reason rather simple models are treated, with a special emphasis on the liquid 
phase. While many of these results have already appeared in one form or another in 
the literature, they are worth repeating. 

The paper is organised as follows: Section || describes the apparent inconsistencies 
that arise between the virial and compressibility routes to thermodynamics for a 
simple density dependent pair potential v(r;p). Section |j] discusses the effective 
pair potentials that result from integrating out three and higher order many-body 
interactions. The effective pair potential that correctly describes the excess internal 
energy is shown to be different from the one that correctly describes the pair structure. 
These points are illustrated with a specific application from polymer solutions. In 
section ^ the McMillan-Mayer|F| tracing out procedure is analysed for an exactly 
solvable lattice version of the Asakura-Oosawajf4| model. While this procedure maps 
onto a useful effective one-component picture in the semi-grand ensemble, integrating 
out the smaller particles in a canonical ensemble does not lead to an effective 
Hamiltonian decomposable as a sum over independent interactions. For charged 
systems the canonical ensemble is the natural choice to integrate out microscopic 
co and countcrions. Again, apparent ambiguities arise when the density dependent 
Debye-Hiickel potential is used to derive thermodynamics. Finally, conclusions from 
these different model calculations are summarised in section |^. 

2. General thermodynamic inconsistencies from a naive application of 
density dependent pair potentials 

As a first rather general example, consider a homogeneous fluid in a volume V, whose 
N particles interact with a spherically symmetric pair potential v(r; p), which depends 
on the global density p = N/V . There are no volume, one-body, or many-body terms. 
No further assumptions as to the origin of the density dependence are made. Two 
established ways to calculate the equation of state (EOS) Z and other thermodynamic 
properties from the correlation functions are[[l5|: 

| Effective potentials could also depend on other state variables like temperature T, but here f mainly 
focus on density dependence 
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(i) the compressibility route: 



(3P f d(]P(p') dp' f dp' 

[l-pc(k = 0;p)\ — , (1) 



P Jo OP P Jq P 
where c(k = 0; p') is the zero-wavelength component of the Fourier Transform (FT) 
of the direct correlation function c(r), [3 = l/fc^T, and P is the pressure. This 
relationship follows from simple properties of the correlation functions and their 
connections to thermodynamics in the grand canonical ensemble - it is therefore 
independent of the particular form of the interactions between the particles, which 
need not be pairwise additive fill, 
(ii) the virial route: 

rP _PP /°° m2 fjv(r;p) Jv(r;p) 



z -=t = 1 --^ p Jo r v^-^^r) 9 ^ (2) 

where g{r) is the radial distribution function. The standard way to derive the virial 
equation is directly through the canonical partition function 

Q(N,V,T) = ^ r J dr N expl-pJ2v( rij -,p)\. (3) 

where A is the usual thermal de Broglie wavelength. The volume derivative in 

^ / 9 iogQ(wn 



dV 



N,T 



also acts directly on the pair potential, which brings in the extra dv(r; p)/ 'dp term in 
the virial equation (||), a result first pointed out in 1969 by Ascarelli and Harrison |fT(| 
in the context of density dependent pair potentials used for modelling liquid metals. 
This particular form of the virial equation is only valid for pair potentials, but the 
derivation of generalisations for systems with three-body terms is straightforward. 

So far so good: both the compressibility equation (which doesn't change from the 
density independent case) and the virial equation (which does) appear to be derived 
for the case of a density dependent pair potential. Nevertheless, this apparent rigour 
deceives, since it is trivial to find density dependent pair potentials where the two 
routes generate different thermodynamics. Consider, for example, a special class of 
density dependent pair potentials with 

v(r;p) = e(p)v°(r). (5) 

Two possible e(p) are shown in Fig. [I]. The compressibility equation (]l|) results in a 
different Z at p = p* for potentials (a) and (b), since the effects of all densities below 
p* are relevant. In contrast, the virial equation (||) cannot distinguish between the two 
potentials at p = p* because it only includes a local density dependence. Of course 
it is not surprising that the two routes to thermodynamics disagree, since one was 
derived in the canonical ensemble, which only samples a single global density, while 
the other was derived in the grand-canonical ensemble, which samples all densities. 

However, the origin of the discrepancy lies deeper than that. One might think 
that the virial equation is less reliable, since it only treats the density dependence as a 
local derivative. But, as the following example will show, the two routes disagree even 
for a simple linear density dependence; even worse, the density dependence correction 
of Ascarelli and Harrison corrects in the wrong direction - it makes the discrepancy 
between the two routes worse. 
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(a) 



Q_ 

to 



(b) 



Figure 1. Amplitude e(p) for the simple potential v(r;p) = e(p)v°(r). 
Potentials (a) and (b) are equivalent except on a small density interval. While 
the compressibility route to the EOS can distinguish between the two potentials 
at p*, the virial route cannot. 



Consider a fluid interacting via the simple Gaussian potential 

v(r;p) = e(p) exp[-r 2 ] (6) 
=e + eip. (7) 

m. 



which could be viewed as a model for interactions between polymer coils 1 17 



This potential falls into the class of mean field fluids ||l9|, |2p| , |2l| , |io|| , for which the 
simple random phase approximation (RPA) c(r; p) « v(r; p) is very accurate and even 
becomes asymptotically exact in the limit of small e(p) or large p. For integrable 
density-independent potentials v(r), the compressibility route leads to an EOS of the 
mean-field (Zmf) form 

Z? PA = Z MF = 1 + ±p/3v(k = 0), (8) 

where v(k) is the FT of the potential. For the same integrable density independent 
potentials, the virial route (|^) reduces to 

Z° mr = Z MF - |#rp jfV {r^|i^} h(r)dr, (9) 

where h(r) = g(r) — 1. In the limit of small v(k — 0) or high densities, the two routes 
to Z approach each other the RPA closure approximation is nearly self-consistent. 

For the density dependent v(r; p) given by Eq. (|^), the compressibility equation (Q) 
takes on a simple form: 

Zr A = l+^(^e p+ 1 -e 1 p^ (10) 

The virial route without density- dependent corrections, i.e. Eq. (||), results in 

Z° mr » Z MF = 1 + tt 3 / 2 Qe p + ieip 2 ) (11) 

for the limit of small v(k = 0) or large p. Within the RPA the corrections to 
Eq. ( pT| ) from the second term in Eq. (|^) can be analytically calculated for Gaussian 
potentials p9| and explicitly shown to be small for the limit being considered. 




Figure 2. Comparison of three routes to thermodynamics for an effective 
potential (3v(r; p) = (1 + O.lp) exp[-r 2 l. Zf PA is given by Eq. @, Z P V ^ PA 
comes from an analytical RPA solution fe 2] to the virial equation fch, and Z® PPA 
from an analytical RPA solution I9J to the simpler virial equation (^|). The 
density-dependent correction to the virial equation increases the disagreement 
between Z c and Z v i r . Inset: In contrast to the density-dependent case, the RPA 
virial and compressibility routes to Z for v(r) = 2exp[— r 2 ] agree within the 
accuracy of the graph, demonstrating the near self-consistency of the closure for 
such potentials. 



At this level, the two routes clearly don't agree. If t\ > then Z® ir > Z c , while 
if ei < then Z° ir < Z c . But, one might argue, the discrepancy should stem from 
ignoring the density derivative term in the virial equation Instead, the opposite 
is true. Adding the density derivative correction results in 

Kir = Z °vir + ^P 2 J r 2 g{rYi exp[-r 2 ]dr. (12) 

Since g(r) > 0, e\ > implies that Z£ ir > Z° ir > Z c , and e% < implies that 
Z% ir < Z° ir < Z c |. In other words, the full virial equation (||), derived explicitly from 
the canonical ensemble with a density dependent potential, is even worse than virial 
expressions (§) or (P) which ignore the density derivative terms. These points are 
illustrated in Fig. |2| for a particular example of the potential (^|). 

The problem clearly lies deeper than the fact that the virial equation only contains 
a local density derivative, since, for the example potential (||), this should be sufficient 
to describe all the density dependence. The next sections provide a partial answer 
to this apparent conundrum by explicitly deriving density dependent potentials, and 
showing that they should really be viewed as mathematical constructs whose physical 
interpretation cannot be separated from the way in which they were derived. 



3. Example 1: effective pair potentials from many-body interactions 

This section demonstrates both in general and with an explicit example that there is 
no unique way to represent averages over many-body interactions as averages over a 

§ The conclusions for the relative ordering of Z v i r and Z c for e\ > 1 hold more generally for all such 
potentials with de(p)/dp > 1 and vice versa. 
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single pair potential that depends on the global density p. Although this has already 
been pointed out several times in the literature]^, |24|, [2^, ^6), mainly in the context 
of noble gasses, it is worth partially repeating because it often seems forgotten. 

3.1. Effective density dependent pair potentials from three-body potentials 

Consider a homogeneous one component fluid system interacting with a density 
independent pair potential (r^, rj) and triplet potential u/ 3 ) (r^, r 3 r^). The 
Hamiltonian can be written as: 

H = K + Y / w {2) (r 1 ,r J )+ £ w& (r,, v 3 , r fc ) (13) 

i<j i<j<k 

with K the kinetic energy operator. By using this Hamiltonian in a standard canonical 
partition function, the following exact expression can be derived for the excess internal 
energy U : 

U{N,V,T) = ^p 2 J dv 1 J dr 2 g 12 w 12 + ^p 3 J dr a J dr 2 J dr 3 g 123 w 123 (14) 

where w 12 = w (2) (ri,r 2 ), w 123 = w (3) (ri, r 2l r 3 ). g 12 = g {2) {r 12 ) and g 123 = 
g( 3 \r\ 2 , 7*i3, r 23) are the homogeneous pair and triplet radial distribution functions 
respectively, and is the distance between particle i and particle j 

Three-body interactions are often cumbersome to use; an effective pair-potential 
which reproduces the properties of the full system would be much more convenient. 
A popular method to achieve this consists of calculating the energy of the system 
governed by the original many-body Hamiltonian through an accurate method (like 
Eq. ( |l4| ) , and then finding an effective pair potential v^f (r^ ; p) that reproduces the 
same energy at the same state-point. For the system at hand, this can be done 
explicitly^, |3J, ||| by writing Eq. Q as 

1 ,2 / . f j„_„ \„,e./ 7 '/ 



U(N,V,T) = -p z J dri / dr 25l2 (r 12 )^ J (r 12 ;p), (15) 

which defines the effective pair potential 

Vu (nj;p) = wi 2 (rij) + Svu(rij',p) (16) 

where 

Svu(rij;p) = \p [ dr 3 g 13 g 23 G 123 w 123 (17) 



is the (additive) density dependent correction to the bare pair potential w\ 2 {r). Here 
G\ 23 is defined as the correction to the Kirkwood superposition approximation ]27| for 
the triplet radial distribution function: 

3123 = 3l2ffl3ff23Gi 2 3. (18) 

In contrast to the virial equation (|^) of the previous section, re-deriving the two- 
body energy equation for a general v(r; p) within the standard canonical ensemble 



route) 15 does not result in extra density dependent terms. In other words, within the 
canonical ensemble, an effective density dependent pair interaction can be constructed 
that correctly captures - by averages over pair correlations alone - the internal energy 
U(N,V,T) of a system with two and three-body interactions. However, Vy (r;p) is 



not sufficient for a complete pairwise description. The three body interaction W\ 23 also 

ef" 



modifies f/i 2 , and although Vy (r; p) generates a different g 12 than would be found by 
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using only W\2 , the pair correlations are not those of the original system, as will be 
shown below. 

A very useful theorem that facilitates the study of pair correlations states that for 
a given homogeneous many-body system at a global density p, there exists a bijective 
one to one mapping between a unique pair potential ifj^ (r; p) and the g^(r) at 
that density |2§|]|j|. While the explicit analytical construction of v^f (r; p) is not as 

straightforward as that of its energy analogue vff^jr; p), an expansion to lowest order 
in density and W123 can be derived |24[ 25| , |3(i| , |3f| : 



fi V12; P) ~ wi2(ri 2 )+Svp(ri 2 ; p) = wi 2 (ri 2 )+p J dr 3 {1 - exp(-wi 23 )} 5i3523-(19) 

Comparing the density-dependent corrections to v e J^ (r; p) and v^j (r; p) in the limit 
of small p and weak W\ 2 ^ leads to: 

To lowest order, the two density dependent corrections differ by a factor of three! 
The unique one-to-one correspondence between g^(r) and Vg^{r;p) therefore 

implies that vff (r;p) cannot reproduce the correct pair correlations for use in 
Eq. (|l5|). This proves that the pair potential derived from the energy equation 
for a many-body system cannot completely and self-consistently capture the excess 
energy within a pairwise description; u| (r; p) is also needed to generate the correct 
g( 2 \r). On the other hand, because the compressibility equation (||) is independent of 
the underlying interactions, Vgf'(r;p) is sufficient to derive the true compressibility, 
and from it other thermodynamic quantities of interest, within a purely pairwise 
description. But the price to pay for coarse-graining a many-body system to an 
effective two-body system in this way is that the thermodynamics can only be 
calculated along one specific route. Neither Vg^ (r; p) nor Vy (r; p) have a well defined 
physical meaning independent of the way in which they were derived. The arguments 
of this subsection are not new, as the following quote, made over 30 years ago in 
exactly the same context, demonstrates: 

We record our opinion that the use of density-dependent effective pair 
potentials can be misleading unless it is recognised that these are mathematical 
constructs to be used in specified equations rather than physical quantities. \2S 



3.2. A worked example: polymers as soft colloids 

The concepts of the previous section can be made more concrete by examining a 
recent coarse-graining of polymers as soft colloids |f7], [l8|, E32j. This system has the 
advantage that many-body correlations and interactions can be accurately calculated 
from computer simulations, allowing detailed comparisons of the coarse-grained system 
with the original many-body system. 

The first step in the coarse-graining procedure is to choose an effective coordinate 
for the polymers, which we take to be the centre of mass (CM). The next step is to 
integrate out the monomeric degrees of freedom to derive an effective interaction 
between the polymer CM. Following the discussion in|Io|, p2[ , the Helmholtz free 

I Note that this pair potential cannot simultaneously represent the correct three-body 
correlations p9| . See e.g. ref. fl32|] for an detailed example and discussion of this point 
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energy T of a set of N polymers of length L in a volume V, with their CM distributed 
according to the set of coordinates {r^}, can be written as the following expansion: 

N N 

F(N,V,{r t }) =^(N,V)+Y,w {2 \r l ,r ] ,)+ £ (r is r„ r fc ) + . . . (21) 

«<j i<j<k 

N 

... + J] IflW^,^...^) 

i<...<N 

In the scaling limit, each term in the series is independent of L as long as the n-tuple 
CM coordinates {ri, r2 . . . r n } are expressed in units of R g , the radius of gyration at 
zero density. This coarse-grained free-energy contains an implicit statistical average 
over all the monomeric degrees of freedom for a fixed set of CM coordinates {r^}. 
J*°\N,V) is the zero-body volume term, related to the internal free energy of a 
single polymer; translational symmetry implies that there is no explicit one-body 
term. The pair and higher body terms are defined in the standard way: the nth body 
term «?(™)(ri,r2 . . . r^J for a particular set of n CM coordinates is given by the total 
free-energy T for n polymers at those coordinates, minus the sum of all the lower order 
termsjl^, |3^] . These interactions can be explicitly calculated by computer simulations; 
examples are shown in Fig. ||, taken from ref. where further details can be found. 
The relative importance of each term decreases for increasing n, so that in principle 
the many-body expansion (^TJ) is expected to converge p2|. 

The free energy F(N, V) of the underlying polymer system follows from a final 
trace over all CM coordinates 

F(N, V) = - log Y, ex P V, {ri})] (22) 

{'«} 

so that Eq. ( |2l| ) can be viewed as an expansion of the effective interaction between the 
CM in terms of (entropic) many-body interactions, in a close analogy to expansion of 
the energy of atomic or molecular systems in many-body interactions f33| . 

But in practise following this path is extremely cumbersome, because the number 
of coordinates and concomitant complexity of the interactions ur") grows rapidly 
with increasing order n. To circumvent this problem one could simply truncate the 
expansion (^) at the pair level, but this would completely ignore the many-body 
interactions. Instead, we recently proposed ]l7|, |T§| , [32| a coarse-graining method which 
takes the many-body interactions into account in an average way. First, at a given p, 
the g{r) between the CM of a polymer was generated by computer simulations. Next, 
for each density, an Ornstein-Zernike integral equation approach was used to invert 
the g(r) and generate the unique [p8| effective pair potential Vg^{r;p) which exactly 
reproduces g(r). Explicit examples of Vg^(r;p) are shown in Fig. ^ for different 
densities p/p*, where p* = 3/(47ri?g) denotes the crossover from the so-called "dilute" 
to the "semi-dilute" regimes|54|. For p — > 0, Vg?f{r\ p) — > w^ 2 \r) while the difference 
5v g (r; p) — v e J ^ (r; p) — (r) grows with increasing density. 

For the polymer system, the expansion ( |l9| ) for the density dependent correction 
to the effective pair potential, Sv g (r; p), can be explicitly calculated ]32|. The results, 
plotted in Fig. [|, show that despite the existence of higher order interactions 
with n > 3 this weak form performs remarkably well, demonstrating that for 
a polymer solution, the three-body interaction is the dominant cause of the density 
dependence in ^■"(r; p) at least for p/p* < 1. 
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Figure 3. Effective potentials «>( n )(r) between the CM of L = 500 self avoiding 
walk polymer coils. The coordinate r denotes the pair distance for n = 2, the 
length of an equilateral triangle for n = 3 and the length of a tetrahedron for 
n = 4 For n = 5 only the interaction at full overlap was calculated, (taken from 
Ref. §§) 
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p/p 
p/p 
p/p 



=0.28 
=0.86 
=2.18 



Figure 5. Explicit simulations of Eq. jis| ) show that the three-body interaction 
(ji, rj i, rfc) is the dominant cause of the density dependence of Vg*(r;p) = 
w^(r) + 5v p (r; p) at lower densities, (taken from ref. fej) 



Since Vg**(r;p) was explicitly constructed to reproduce the correct pair 
correlations, it follows that the true thermodynamics of the full many-body system 
should be reproduced by using this potential in the compressibility equation (|l]). The 
results shown in Fig. (JsJ) confirm this: full polymer simulations of the EOS Z for 
L = 500 and L = 2000 self avoiding walk polymers on a cubic lattice compare well with 
the EOS Z c calculated with the appropriate (r; p). The small residual differences 
at large p are most likely due to numerical difficulties in performing accurate inversions 
at these high densities |J5|. 

We now have a concrete example where the explicit construction of Vg**(r;p) 
leads to the full thermodynamics of a many-body system within a purely pairwise 
description. But what happens when this density dependent pair potential is used 
within the virial equation (^)? Fig. |^ shows that the virial equation without 
the explicit density derivative is accurately approximated by the mean-field form 
Z° ir w Zmf = 1 + pv(k = 0;p). Just as was found for positive t\ in section Z® ir 
overestimates the true Z, and adding the density-dependent correction (not shown) 
results in a significantly larger overestimate. 

These results provide a partial explanation for the thermodynamic inconsistencies 
found in section g, If v(r; p) is equivalent to v e g ' fJ r (r; p) then only the compressibility 
equation has an unambiguous physical interpretation. On the other hand, one can also 
follow a route similar to that used to derive v^/^(r; p) to obtain a p) which will 

reproduce the true pressure of an underlying many-body system through the virial 
equation g6[ |36j] . This potential will not equal Vg^(r;p), and so won't generate 
the correct pair correlations needed in Eqs. (||) or (^). For an underlying many-body 
system, no single pairwise density dependent pair potential v(r; p) exists which can, 
through either the full density dependent (||) or the density-independent (Q) forms of 
the virial equation, yield the true thermodynamics. 

Nevertheless, Fig. ^| shows that the relative overestimate found using Vgff(r; p) 
in the virial expression is less than the relative underestimate found when all many- 
body terms are ignored by only taking w/ 2 )(r) into account. Using u|^(r; p) in the 
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Figure 6. Equation of state of a polymer solution calculated in several ways: 

Circular symbols are_fpr direct simulations self avoiding walk polymers; the two 

Z c follow from Eq. (II); within the statistical errors of the method, they agree 

with the direct simu ations of the full polymer system. Z°. was calculated by 

' J J vvr J 

computer simulations with Vg^(r; p) Jlq| , and agrees well with the simpler Zmf- 

In contrast to Z c , this is not a good estimate for the true EOS. Finally, when 

Zmf is calculated without any density dependence of the pair potential, i.e. with 

«j( 2 )(r) only, it strongly underestimates the EOS. 



simple virial expression ( pi] ) is therefore a better approximation than either ignoring 
density dependence all together (as with w^ 2 \r))^ or using the full density dependent 
form (|2|) of the virial equation. In fact for the dilute regime, p/p* < 1, where 
the polymers as soft colloids coarse-graining technique is most useful, the absolute 
differences between the virial and compressibility routes to the EOS are quite small. 
Differences in the structure generated by the two effective potentials are also rather 



small) 35 . On the other had, for the semi-dilute regime (p/p* > 1), the polymer 
EOS scales as Z cx p 1 /^-!) ss p 1 - 3 , where the Flory exponent v w 0.59@. At 
these higher densities the RPA approximation is excellent, which suggests that the 
dominant density dependence of the FT of Vgf*(r;p) scales as u|"(0;p) ~ ap 3 , 
with a a density independent constant. This implies that Z c ~ (a/2.3)p , while 
Z® ir w Zmf ~ (a/2)/} 1 " 3 . Using v|"(r;p) in the virial equation at high densities 
results in the correct scaling exponent, but a prefactor which is about 15% too high. 
In contrast, ignoring the many-body terms altogether by using only io' 2 )(r) in either 
the virial or the compressibility equations, results in Z ~ p, which scales with the 
wrong exponent. 



3.3. Lessons for coarse- graining many-body systems 

The two previous subsections have shown that there is no unique way to represent the 
effects of many-body interactions in effective pair interactions. This has important 
implications for several techniques to derive such effective pair potentials. For example, 
if an ab-initio or other higher level approach is used to generate the structure of a 
many-body system, then the effective pair potential that reproduces this structure will 
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not be the same as the pair potential that correctly describes the internal energy. When 
the many-body interactions are weak, these differences may not be that important, 
but when they are strong, they may be significant. However, as seen for the polymers 
as soft colloids approach, using Vg^(r;p) in the virial equation is still a better 
approximation than completely neglecting all many-body interactions. One might 
hope that the same is true for methods that derive effective pair potentials from 
internal energies or from structure. 

In general, the stronger the many-body interactions, the more likely that extra 
care must be exercised when applying an effective pair potential derived by one 
route (i.e. structure) to extract other physical quantities (i.e. energy). This has 
implications for work done on systems with strong angular forces like water, where 
rather complicated effective pair-potentials have be constructed to mimic certain 
properties gj. 

A more physically motivated way to average over many-body interactions may 
be to derive effective pair interactions with a local density dependence, since when 
the local instantaneous density of a liquid is higher than the global average, one 
expects the relative strength of the many-body interactions to be more important 
there, and vice versa for lower local densities. In fact, for a number of systems, 
one can show explicitly that a local density dependence is equivalent to many- 
body interactions |33|, Q (the same can be done for internal degrees of freedom Q). 
Although this doesn't imply that any many-body interaction can be consistently 
mapped onto a local density dependence, for those cases where it can be done, 
a completely self-consistent thermodynamics should exist based on these effective 
potentials. 

Another obvious problem with a global density dependence arises when one tries 
to treat inhomogeneous systems. Even for Argon, the LJ pair potential generates a 
surface tension that differs by up to to 19% compared to calculations that include 
explicit three-body effects p2|. Again, a local density dependence would appear the 
more natural way to coarse-grain. Of course this opens up new problems, such as how 
exactly does one define the local density etc. . .fij 



4. Example 2: effective one-component Hamiltonians for two-component 

systems 

Effective state dependent potentials naturally arise from tracing out one component in 
a two-component mixture. For uncharged systems such procedures were first carefully 
formulated by McMillan and Mayer in their famous theory of solutions|ll|. As they, 
and many other authors have stressed (see e.g. j2?| |]| Q), this works most naturally 
in an ensemble where the component to be traced out is treated grand-canonically. 
For charged systems this is no longer the case; instead, charge-neutrality makes the 
canonical ensemble the natural choice. 

As the following pedagogical examples will illustrate, whereas tracing out 
procedures for uncharged systems are quite well understood, for charged systems the 
waters are still muddied. 

% An interesting recent example for an NVT simulation can be found in^oj. However, the original 
effective potential studied was of the v y ff (r;p) formQ, while in K(t it is used in a virial equation 
like Eq. (fcj) . The lack of thermodynamic self-consistency found by these authors stems in part from 
their use of a density dependent interaction without a careful analysis of its origin. 
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Figure 7. A 2-dimensional lattice version of the AO model. The big squares 
live on the N sites, and the small ones live on the 2N links. The big squares 
are hard, and cannot overlap each other, and similarly the big and small squares 
exclude each other. The small squares are ideal w.r.t. each other, and so multiple 
occupancy of links not excluded by big particles is possible. The crosses on the 
central big squares show which links they exclude from the small particles. 



4-1- Asakura Oosawa and related models 

A popular model for polymer-colloid mixtures consists of treating the colloids as simple 
hard spheres (HS) of radius R c and the polymers as penetrable sphere (PS) whose 
interaction with the colloids is HS like, with a cross diameter o~ cp of order R c + R gi 
but whose interaction with other the polymers is ideal gas like JT^| . Here I introduce a 
very simple and exactly solvable lattice version of this Asakura Oosawa (AO) model, 
described in Fig. [?]. 

In the grand canonical ensemble the partition function for the mixture becomes 

E m (23) 

{rii} {nii} \ i i / 

where [ii denotes the chemical potential of each species, and the sums are over all 
possible configurations {ni} and {mj of the big ("colloidal") and small ("polymeric") 
particles respectively. Since the model is athermal, one can set [3 = 1; in addition, 
the de Broglie wavelength is set to A = 1 which simplifies notation throughout this 
section. 

For a given configuration of the big particles {ni}, there will be M({ni}) links 
left for the small particles. Their partition function can then be calculated 

M M l z l 

H S ma«({M) = ^ —fT = exp[z s M({n. i })] (24) 

1=0 

where the fugacity z s = exp[/3/i]. For a fixed configuration {rii} of big particles, the 
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number of free links is: 

M({ni}) = 2N - 4 rk + E Ul n i ( 25 ) 

i <ij> 

since there are two links for each site, each big particle excludes 4 links, but when 
two big particles touch, they only exclude 7 links, because one is doubly excluded, as 
illustrated by the crosses in Fig.|. This increase in the effective volume available (and 
concomitant entropy) for the small particles when two big ones touch is the origin of 
the osmotic depletion effect ]l4| . A large enough fugacity of the small particles can 
induce an entropy driven phase-separation of the big particles. 

With the explicit form (^5|) for the available "free volume" M{{ni\), the full 
partition function (|23|) simplifies to 



Zmix(Hb, fi s ,N) = exp [2z s N] Y ex P 

{«*} 



(/3/Z 6 - 4z s ) Y n i+ Z s Y Ui n 3 



<ij> 



.(26) 



This illustrates the essence of the McMillan-Mayer tracing out procedure: the sums 
over states in the partition function of the original two-component system are rewritten 
as a single sum over the states of an effective one-component system. There is no need 
to invoke any Born-Oppenheimer type time-scale separation arguments; the procedure 
might even be viewed as simply a mathematical trick. One could just as well trace 
out the big particles. But the particular advantage of tracing out the small particles 
is that the effective one-component system is simpler to treat, since it can be mapped 
onto a binary lattice gas model with an effective chemical potential /z^ — fib — 4z s 
and an effective Hamiltonian that, in turn, is decomposable as a pairwise sum over 
nearest neighbours with an interaction strength (3e e ^ = —z s . By mapping onto an 
Ising spin system in 2-D, this model, with its entropically driven phase-separation 
transition, can be exactly solved, as shown previously for a closely related model of 
non- additive hard squares^]. Mathematically these models are a special cases of the 
more general decorated Ising model introduced by Widom[[f6"|. 

Keeping both species grand-canonical facilitates the mapping to an Ising model, 
but it is just as easy to integrate out the small particles in a semi-grand ensemble, i.e. 
fixing N b , z s , N, resulting in an effective partition function of the form: 



,(N b , z s ,N) = exp [z s (2N - AN b )} ^ exp 



E < 

<ij> 



(27) 



where {n^}' denotes all possible ways of arranging N b particles on a lattice of N sites. 
Again, as long as z s is fixed, Z m i X can be interpreted as an effective one-component 
system, interacting with an effective Hamiltonian of the form 

pH e "(N b: z s ,N; {n,}) = H bb - z s 2N + z s 4N b - z s £ n t nj (28) 

<ij> 

with a configuration independent "volume term" , and a pairwise decomposable pair 
interaction. H bb is the bare hard-core big-big interaction. All the standard statistical 
mechanics for such one-component systems can now be brought to bear to calculate 
correlations and phase-behaviour. The volume term contributes to thermodynamic 
quantities, but not to the phase-behaviour fill M. 
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The tracing out procedure can also be done in the canonical ensemble, keeping 
both N b and N s fixed: 



Z mix {N b ,N s ,N) 



{m}' 



<ij> 



(29) 



{m}' 

But now the effective one-component system is no longer equivalent to a system with 
a pair-decomposable effective Hamiltonian; at best it can be rewritten as: 



H e " (N b , N s , N; {n t }) = H bb - N s log [2N - 4N b ] - N s log 



1 



1 



2N - AN, 



(30) 



(31) 
(32) 



In the limit of only two big particles the logarithms can be expanded and pair term in 
the Hamiltonian reduces to the same form as found in the same limit for the grand- 
canonical tracing out procedure. But for a larger number of big particles this ceases 
to be true; the Hamiltonian can no longer be written as a sum over independent many- 
body interactions of the form of Eqs. @ or @. The effect on H ef f (N b , N s , N; {n-}) 
of changing the number of pairs in a configuration {n{\ by one depends on the 
configuration of all other pairs in the system. The McMillan-Mayer mathematical 
tracing out procedure does not lead to such a useful simplification in the canonical 
ensemble as it does in the semi-grand ensemble. 

Similar manipulations can be performed for off-lattice two-component systems 
such as the original AO model Jl4|]. Integrating out the N p PS polymeric particles for 
a fixed configuration jr,} of N c colloidal HS results in effective Hamiltonians of the 
form: 

H eff (N c , z p , V- {r, ; }) = H cc + Q(N C , z p , V- {r,}) 
H e "{N c , N p , V; {r, ; }) = H cc + F{N C , N p , V; {r J) 

for the semi grand and canonical ensembles respectively. H cc is the bare colloid 
Hamiltonian. £l(N c , z p , V; {r^}) is the grand-potential and F(N C , N p , V; {r^}) is the 
Helmholtz free energy of an inhomogeneous system of PS particles in the external 
field of the colloids. The free energy form of these effective Hamiltonians holds for a 
general two-component system[|4). Integrating over all possible configurations {r^} of 
the colloidal particles leads to the full two-component partition function. For the AO 
model these Hamiltonians take a particularly simple form since 

n(N c , z p , V; {r,}) = -z p V free (N c , V; {r,}) (33) 

is the grand potential of an ideal gas at in the accessible free volume V^ ree (N c , V; {r,}) 
and similarly 

F(N C , N p , V; {r,}) = N p \og[N p ] - N p - N p log [V free (N c , V; {r,})] • (34) 

is the Helmholtz free energy of an ideal gas of N p particles in the accessible free volume. 

The calculation of V^ ree (N c , V; {r^}), the direct analogue of M({ni\) in the AO 
lattice models, simplifies dramatically if the size-ratio R g /R c < 0.1547. There are 
then no triplet or higher order overlaps of the exclusion zones, and the accessible free 
volume can be written as 

V^{N C , V; {ri}) =V-N c Vx+J2 ^a(r«) 

i<j 



(35) 
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where V\ = |7r(cr cp ) 3 is the volume excluded by each colloidal particle, and V 2 (r) has 
the standard AO formal), depending only on the relative separation r of two particles. 
(Note the correspondence with Eq. j25| ) for the lattice model, where V = 2N, Vi = 4 
and the sum over V2M is replaced by a nearest neighbour lattice sum). The effective 
Hamiltonians of Eqs. (pll) and (p2|) simplify to 



H e ff(N c , z p , V; {r 4 }) = H cc - z p V + z p N c V x - z p V V 2 {r) (36) 



> / 

i<j 



H e "(N c , N p , V; {r t }) = H cc + N p \og\ %— ] - N p 



L V-N c Vi- 



N p log 



V - N C V! ^ 



2 P V 2 (r) 



(37) 



Just as was found for the lattice model, the canonical tracing out procedure leads to 
a Hamiltonian that cannot be written as a sum over independent interactions of the 
form of Eqs. © or @)Q In contrast, semi-grand H e ^ (N c , z p , V, {r^}) takes on the 
more useful form of a "volume term" plus a sum over an effective pair potential. The 
interpretation of the latter will now be treated in more detail. 

While volume terms in Eq (^) are important to make contact with the full two- 
component system, as worked out very clearly and in great detail by Dijkstra van 
Roij and Evans E^J, they do not contribute to the phase-behaviour of the system. For 
Rg/R c < 0.1547 the simplified pairwise Hamiltonian 

H e ^{N c , z p , V) =H CC + J2 < ff {r\ z P ) (38) 

i<j 

with v e Jf(r\z p ) — —z p V2(r), completely determines the phase-behaviour and 
correlations of the colloidal system. An entirely self-consistent one-component osmotic 



thermodynamics can be defined based on the Hamiltonian (38). This means, for 
example, that the compressibility equation ([j]) generates the same osmotic EOS as 
the virial equation: 



Pc 3 J Or 

where n(A^ c , z p , V) is the osmotic pressure^. Since z p is a fixed external parameter, 
v e J^(r;z p ) is in fact not state-dependent within this one-component picture. Thus, 
the McMillan-Mayer tracing out procedure results in a system that can be exactly 
mapped onto a classical one-component fluid with a state-independent pair potential]. 
One must simply keep in mind that the resultant thermodynamics quantities are 
osmotic, and defined w.r.t. a reservoir containing PS particles only. Similarly, since the 
McMillan-Mayer procedure is essentially a mathematical trick to facilitate calculating 
the sums in the partition function, very little information can a-priori be extracted 
about the effective dynamics of the one-component system. For that, further physically 
motivated arguments and approximations are needed. 

+ The small differences between the volume terms in the lattice and off-lattice versions of the AO 
model stem from kinetic energy terms that are ignored in the former case. 

f Note that adding the explicit volume term dFo/dV to the compressibility and virial equations 
would generate the thermodynamics of the full two-component system. 

I For larger size-ratios, the mapping results in three-body and higher order interactions, leading to 
similar problems to those discussed in section H when one tries to derive an effective pair potential. 
Also, interactions between the small particles will lead to effective three-body and higher terms. 
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Now suppose one were working in the canonical ensemble, and naively used the 
Asakura-Oosawa depletion potential (r; z p ), as has sometimes been done in the 
literature. Then z p = z p {p c , p p ) would no longer be an external parameter but would 
instead depend on the state of the system. For example, for small size-ratios a good 
approximation is z p sa /0 P /(1 — /9 c ^i)@, ffTj- One might be tempted to define an 
effective one-component thermodynamics by fixing p p , and treating it as an external 
parameter, as might be natural in an experiment. The potential then takes on the p c 
dependent form v e Jf (r; p c ; p p ). If this potential is taken as given, without enquiring 
as to its origins, as was done in section |2|, then the compressibility and virial equations 
take forms similar to Eqs. (0) and (g) respectively, since v^f? depends on p c . One might 
expect that the ensuing thermodynamics would again be that of the osmotic system, 
but now neither the derived virial nor the derived compressibility equation is correct. 
For the present case of a fixed p p , the compressibility equation (|l|) using v p (r; p p , p c ) 
does not result in the correct osmotic pressure because z p varies when p c is changed. 
On the other hand, using v e p ^ (r; p c , p p ) in the simpler virial equation (|39|), i.e. one 
without a density derivative, does generate the correct osmotic EOS, as follows from 
the following arguments: The thermodynamic properties are independent of ensemble, 
so, for a given state point (N c , N p , V), the osmotic pressure is the same as in the semi- 
grand ensemble at a state point (N c , z p , N) such that < N p > zp .Nc,v= N p . There the 
potential vyf (r; p p , p c ) = v^^(r;z p ) generates the correct pair correlations g cc (r), 
and also the correct virial pressure U(N C , N p , V) through Eq. (f^). But the apparent 
relevance of this osmotic virial equation within the canonical ensemble is deceptive - 
it only follows because it can be derived in the semi-grand ensemble, and used at a 
state point where v e J' (r; p c , p p ) = u|" (r; z p ). 

Of course the lack of consistency between virial and compressibility routes should 
not be surprising, since a careful tracing out procedure demonstrates that the effective 



Hamiltonian (37) for the canonical ensemble cannot be decomposed into a sum over 
independent pair potentials. Therefore, using the Asakura-Oosawa pair potential in 
this ensemble is not rigorously justified, except for instances where parallels with the 
semi-grand ensemble can be madc|| 

In summary then: for a fixed z p , a completely self-consistent one-component 
thermodynamics can be derived in the semi-grand ensemble for an AO system. But 
within the canonical ensemble, a McMillan-Mayer style tracing out procedure leads to 
an effective Hamiltonian that cannot easily be written as a sum over pair and higher 
order interaction terms; the ensuing one-component system does not have a simple 
interpretation as an effective liquid. Moreover, if the AO depletion potential is naively 
applied in the canonical ensemble for fixed p p , then the density-dependence of the 
pair potential again leads to an apparent lack of consistency between different routes 
to osmotic thermodynamics. In contrast to the case studied in section ||, where the 
effective density dependent potential Vg^(r; p) that arises from tracing out three-body 
interactions generates the correct thermodynamics only through the compressibility 



equation, here only the virial equation (|39|) should be used 



§ At this point I should point out that in the limit N c = 2, both the canonical and the grand- 
canonical ensembles do result in the same pair potential, as can be seen by expanding the log in 
Eq. (|37|). A similar conclusion holds for higher order interactions, should they be relevant. However, 
I stress that great care must be used when these low density effective interactions are applied at finite 
colloid densities in the canonical ensemble. 
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4-2. Debye Hiickel model 

The effective interactions and resulting phase-behaviour of charge-stabilised colloidal 
suspensions have been the subject of much recent debate^, In contrast to 
uncharged mixtures, global charge-neutrality implies that the canonical ensemble is 
the natural choice in which to integrate out the co and counter ions, to arrive at an 
effective one-component colloidal picture. And this explains in part why the problem is 
so difficult, since, as was shown in the previous subsection, tracing out one component 
in the canonical ensemble does not necessarily lead to an obvious description in terms 
of independent (many-body) interactions. In addition, direct computer simulations 
of the full mixture are greatly complicated by the long-range nature of the Coulomb 
interactions and the large length scale differences between a typical colloidal particle, 
and the co- and counter-ions. 

Rather than attempting yet another tracing out procedure, this section has a 
much more modest goal, namely to illustrate pitfalls that arise from a naive application 
of a very simple textbook density dependent potential of the Debye-Huckel screened 
Yukawa form 

Z 2 

l3v DH (r; p) = — exp [-n(p)r] . (40) 
r 



Here Z is the charge of the colloidal particle, and n(p) = y/AirZp is the screening 
parameter in the absence of salt. The Bjerrum length Xb ~ f3e 2 /e, with e the elemental 
charge and e the dielectric constant, has been set to l,to simplify the notation. Since 
k depends on the overall density (through charge neutrality), it should come as no 
surprise that a simple application of the compressibility equatio n p ) and the virial 
equation (gj) do not generate the same thermodynamics. Since Eq. (|40| ) is an integrable 
potential, 

Pv(k = 0;p) = ^ = -, (41) 

K p 

its thermodynamic behaviour resembles that of a mean field fluid (h]] for large p or 
small effective Z, where the RPA closure should be quite accurate. Thus the two 
routes lead to: 

z rpa = 1 + z> (42) 
which can be interpreted as the ideal EOS of the colloids and Z counterions, and 



Z vir = (1 + ~Z) - ~/?7rp J t 2 



, , w dv DH (r) dv DH (r) 
h ( r )( r aZ ) -fW3p- 



(43) 



dr dp 

Even the leading term in the virial equation differs from the compressibility equation. 
Since both the r and p derivatives of VDH(r;p) are always negative, the second two 
terms of Z V i r both reduce its value w.r.t. the leading Zmf = i^ + ^Z) term|], increasing 
the difference between the two routes even further. The present discrepancy originates 
not in the lack of consistency of the closure, but rather in the naive application of a 
density dependent pair potential. A more careful analysis of the underlying two- 
component colloid + counterion system shows that volume terms must also be taken 
into account p8| , but these don't bring the two routes any closer together. The only 

| At least in the regime where the RPA is a reas ona ble approximation. Note that the virial equation 
can be exactly solved in the RPA approximation |22| 
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way to know which (if any) of the two routes is the more reliable would be to derive 
them from a careful analysis of the full tracing out procedure^ 

The arguments above may have repercussions for the modelling of liquid metals by 
effective density-dependent pair potentials. There, the Ascarelli-Harrison form of the 
virial equation (||) has often been applied, but as other examples in this paper show, 
it is not yet clear whether this equation is reliable. In fact, a recent careful derivation 
of an effective one-component virial equation within within a two-component electron 
ion picture by Chihara et a/.p9]| does not include an explicit density derivative. 

Finally, it appears that effective potentials in charged systems will always be 
density dependent due to charge-neutrality. But this same charge neutrality implies 
that a two component system can be viewed effectively as a one-component one. For 
example, all the partial structure factors have the same k — > limit, up to trivial 
prefactors, i.e. they are not independent. The density dependence in charged systems 
most likely has a different character from that found in uncharged systems. 

5. Conclusions 

The overall conclusion of this paper is simply that 

an effective density dependent pair potential v(r; p) cannot be 
properly interpreted without reference to the coarse-graining 
procedure by which it was derived. 

This was shown by a number of explicit examples. The only systems where an effective 
pair potential could be rigorously interpreted as part of an effective Hamiltonian were 
the AO model for R g /R c < 0.1547 and its lattice analogue. Other potentials like 
Vgff(r; p) cannot be formally interpreted as part of an effective Hamiltonian. Instead, 
they should be viewed as mathematical devices to calculate desired properties from a 
particular equation |23|, pL |2f|. 

Some of the more detailed points are repeated below: 

• A pair potential that depends on global density, v(r; p), does not lead to the same 
thermodynamics in the canonical and grand-canonical ensembles. 

• The status of the Ascarelli-Harrisonfl^l form of the virial equation (Q), with 
an added dv(r; p)dp term, is suspect. For uncharged systems and many- 
body systems, counter examples where it does not apply and even makes 
thermodynamic inconsistencies worse, can easily be constructed. Whether it is 
valid for liquid metals or other charged systems remains to be proved. 

• There is no unique way to represent the effect of many-body interactions as 
density dependent pair interactions (23[ |25|, p6| . The potential that correctly 
reproduces the structure, vy'(r;p), will not generate the right internal energy 
when used in the energy equation etc... At best, one can pick a v e ^(r;p) that 
performs well for the particular physical properties one is interested in. For 
example, the parameters commonly used in the LJ potential to model Argon can 
be viewed as a compromise between those that correctly reproduce the energy, 
the virial pressure, and the structure p^]. 

One example where Eq. (^) appears to be better than forms without the density derivative can be 
inferred from the discussion of the Debye-Hiickel model by Belloni[p|. By adding the explicit volume 
term to his application of the Debye-Hiickel potential in Eq. (^), the virial pressure closely resembles 
an approximate form derived from a two-component picture. 
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• The previous statement implies a link between problems with the transferability of 
an effective pair potential, and problems with its representability. For example, if 
many-body interactions generate a large relative density dependence in v e ^ (r; p), 
then the potential derived by a given route at p\ will differ significantly from the 
one derived at pi (transferability problems). This large density dependence also 
implies important differences between v^j (r; p) and Vgff(r;p) (representability 
problems) . 

• A McMillan-Mayer type tracing out procedure for a two-component system in 
the canonical ensemble docs not lead to an effective one-component Hamiltonian 
that be written as a sum of independent interactions. 

• The Debye-Hiickel potential, like all such density dependent pair potentials, leads 
to thermodynamic inconsistencies when it is naively used in an effective one- 
component picture. At present it is not yet clear what the correct route to 
thermodynamics is. 

• Obvious problems arise when potentials like v(r; p) are used in inhomogeneous 
systems. There it would make more sense to coarse-grain to potentials that 
depend on a measure of the local density. In some cases, this can be shown to be 
equivalent to a many-body interaction approach. 

I finish with the question: when do these results matter? The picture is 
not as bleak as I might have painted it in the discussions above. Indeed, in the 
case of Argon, the differences between Vg^(r;p), and v^^{r;p) are not so large, 
so that a single compromise effective pair potentials works admirably well for the 
liquid phase. In addition, many examples can be found for soft-matter systems 
where ignoring the many-body forces altogether, - simply using w^ 2 \r) in the 
expansion of an H e ff - works quite well. This is true for many systems described 
by depletion potentials |5(], [44], pl| , or for the effective interactions derived for star- 
polymers Q or dcndrimcrs|53|. And in other cases, such as the polymers as soft 
colloids example of section [|using vt^(r;p) in the simple virial equation may still 
be a better approximation than ignoring many-body effects altogether. Nevertheless, 
the particular coarse-graining procedure used to derive v e ^ (r; p) must always be kept 
in mind - the naive consumer of effective pair potentials should beware. 
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